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^  An  Integrsl  represents! ion  technique  is  presented  for  determining  the  mods! 
propagation  properties  of  a  homogeneous  cylindrical  dielectric  waveguide  of 
arbitrary  cross-sectional  shape  and  Indes  n  j  embedded  In  a  medium  of  Index 
n_.  Result  a  are  given  for  weakly  guiding  fibers  of  various  shapea.  Among 
tlfese  are  re<-**ngles  and  ellipses,  which  makes  comparisons  with  previous 
work  possible,  w 


tnclasalfled 


Summary 

A  technique  is  presented  for  determ ining  the  modal  propagst i<«n  properties  of 
a  homogeneous  cylindrical  dielectric  waveguide  of  arbitrary  croaa -sectional  shape 
and  index  nt  embedded  in  a  medium  of  index  n^.  For  the  cylinder  parallel  to  the 
r-axta  atl  field  components  can  be  derived  from  E^  and  bf.  Integral  repreaenta- 
tlona  for  E^  and  By  are  derived  which  aattafy  the  appropriate  Helmholt7  equations 
inside  and  outside  the  guide  and  which  guarantee  that  the  boundary  conditions  are 
satisfied.  On  expanding  Ef  and  B?  in  certain  sets  of  basis  functions  the  integral 
representations  become  a  set  of  linear  equations.  The  vanishing  of  the  determi¬ 
nant  of  this  set  yields  the  propagation  Constanta  of  the  various  modes.  For  the 
important  special  case  of  weakly  guiding  Tiber*  (nt  •  n2>.  Ey  and  B^  become  small 
and  we  deal  instead  with  the  relatively  large  transverse  rectangular  components 
of  E  and  B.  any  one  of  which  satisfies  a  Helmholtr  equation  inside  and  outside  the 
guide.  A*  with  the  general  rase,  an  Integral  representation  is  derived,  basis 
functions  are  inserted,  and  a  determlnantal  equation  is  generated  whose  roots 
yield  the  modal  propagation  constants.  Iteaults  are  given  for  weakly  guiding  fibers 
of  various  shapes.  Among  these  are  rectangles  and  ellipses,  which  makes  com¬ 
parisons  with  previous  work  possible. 
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Mod**  of  Dioloctric  WoveguicUs 
of  Arbitrary  Cross  *  S*ctional  Shop* 


i  i  vtikhh  cno* 

With  the  increased  interest  tn  fiber  optic  or  integrated  optic  communication 
technique*  haa  come  a  need  to  know  in  tncreaatng  detail  the  propagation  propertiea 
of  cylindrical  dielectric  waveguide*  (Gloge,  *  Goell,  Marcatlll3).  In  fact,  the 
literature  haa  so  burgeoned  that  it  would  be  excessive  to  try  to  quote  all  the 
individual  paper*  that  discuaa  the  subject.  Happily,  there  are  by  now  aeveral 
excellent  books  and  review  article*  that  make  this  unnecessary  (Marcuse,  * 

Kapanv  and  Burke5). 

The  problem  of  electromagnetic  propagation  down  a  dielectric  cylinder  is 
closely  related  to  the  problem  of  scattering  of  electromagnetic  waves  from  the 
same  cylinder.  It  is  then  not  surprising  that  methods  used  for  the  scattering 
problem  <Eygesfi>  have  also  been  applied  to  the  problem  of  propagation.  Thus, 

(Received  for  publication  3  July  1979) 

1.  Gloge.  D.  (1971)  Appl.  Opt.  10:2252. 

2.  Goell.  J.E.  (1969)  Bell  Svs.  Tech.  J.  48:2133. 

3.  Marcatili,  E.A.J.  (I960)  Bell  Sys.  Tech.  J.  48:2071. 

4.  Marcuse,  D.  (1974)  Theory  of  Dielectric  Optical  Waveguides.  Academic 

Press.  N.  V. 

5.  Kapany,  N.S.  and  Burke,  J.J.  (1972)  Optical  Waveguides,  Academic  Press. 

N.  Y. 

6.  Eyges,  1..  (1973)  Ann.  Phys.  81:567. 
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point  matching  (Goell  )  and  ray  tracing  (Love  and  Snyder;  de  Vita  )  technique* 

have  been  applied.  Deaplte  much  work,  however,  there  la  no  one  general  analytic 

method  which  appltea  to  atngle  guide*  of  arbitrary  croaa-aect tonal  ahape  and  to 

the  coupling  between  two  or  more  auch  guide*. 

In  thia  paper  we  preaent  auch  a  method  for  aingle  guide*;  ita  extension  to 

coupled  guide*  la  in  preparation.  The  method  la  largely  analytic  although  finally 

machine  computation  la  requir  d.  The  method  not  only  provide*  new  technique* 

for  irregularly  ahaped  guide*,  but  alao  haa  the  minor  advantage  of  providing  neat 

derivation*  of  the  standard  formulae  for  atngle  circular  guides,  including  the  well- 

known  approximation*  that  ariae  when  the  index  of  refraction  i*  little  different 

from  that  of  the  surrounding  medium.  The  preaent  method  is  an  extenaion  of 

recent  new  technique*  for  solving  the  problem  of  scattering  by  irregularly  ahaped 

dielectric  bodies,  and  in  the  static  limit,  for  aolving  the  problem  of  an  irregular 

6  8 

dielectric  or  permeable  body  in  an  external  field  (Eyges,  Eyges  and  Gtantno, 

9 

Nelson  and  Eyges  ).  The  general  idea  here  la  the  *atne  aa  in  those  paper*.  That 
is,  an  integral  repre*entat ion  for  the  various  field  component*  i*  derived  that  ta 
equivalent  to  the  appropriate  differential  equation*  in*ide  and  outaide  the  guide  and 
that  guarantees  a*  well  that  the  various  boundary  condition*  are  satisfied.  These 
differential  equation*  for  the  field  components  inside  the  guide  have  a  set  of  known 
solutions  that  can  be  considered  aa  a  set  of  baaia  function*.  The  field  component* 
can  be  expanded  in  a  aerie*  of  these  basis  function*,  with  unknown  coefficient*, 
on  introducing  the  expansions  into  the  integral  representation*,  they  become  a  set 
of  linear  (matrix)  equations  whose  matrix  element*  involve  theae  basis  function* 
in  line  integral*  that  are  taken  over  the  boundary  of  the  cross-sectional  shape. 

This  is  the  only  way  the  boundary  enters;  there  is  no  need  to  match  interior  and 
exterior  solutions  arms*  it.  Since  the  set  of  matrix  equations  Is  a  homogeneous 
one.  it  ha*  solutions  only  for  certain  allowed  value*  of  the  parameters  Involved, 
t *ie  of  these  parameter*  is  the  propagation  constant  we  seek. 

The  description  above  applie*  to  a  guide  of  index  nj  embedded  in  a  medium  of 
smaller  index  n^  but  with  no  other  restrictions  on  the  magnitudes  of  nj  and  nj.  In 
practice  however,  it  is  very  common  for  nj  to  be  close  to  n^.  In  this  "weakly 
guiding"  rase  rertain  vectorial  aspects  of  the  general  procedure  outlined  above 
become  simpler,  and.  in  fart,  the  problem  ran  be  reduced  in  a  good  approximation 
to  a  scalar  one.  We  have  therefore  concentrated  on  thia  scalar  formulation  and 
have  based  our  numerical  calculations  on  it. 

7.  Love.  J.D.  and  Snyder.  A.W.  (1P77)  Ann.  Telecommunications  32:100; 

de  Vita,  P.  (1077)  Ann.  Telecommunications  32:1  Iftl 

8.  Eyges,  L.  and  Giant  no.  P.  D.  (1079)  IEEE  Tran*.  Ant.  <L  Prop.  AP-27:5!)7. 

9.  Nelson.  A.  E.  and  Eyges,  L.  (1078)  J,  Opt.  Soc,  Am.  88:254. 
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2.  INTEGRAL  REPRESENTATIONS  K)K  THU  NEARLY  GUIDING  CASK 


In  thta  section  integral  repreaentationa  in  the  weakly  guiding  caae  are  derived 

for  the  cylmdrica*  guide  ot  arbitrary  rro»«  aection  aketched  in  Figure  1.  The 

time  dependence  ta  given  aa  e*u*’,(  ao  that,  for  example,  the  electric  field  E(r,  t) 

ta  written  aa  E  E  e'w’’*.  The  guide  medium  la  as  Burned  uniform  with  index  of 
—  —  w  2 

refraction  nj  and  dielectric  conatant  *  j  related  by  nj  •  •  j  »nd  aimilarly  the 

external  medium  ia  characterized  by  02  «  2-  Then  E  (aimilarly  H  )  aatiaflea 

< V*  ♦  k.)  E  0  inaide  the  guide  and  (rz  ♦  kf)  E  •  0  in  the  external  medium. 

where  k,  '  n.k  ,  k,  n,k  and  k  •  w/c. 

1  1  o*  2  2  o  o 


y 


Figure  1.  Croaa  Section  of  a  Cylindri¬ 
cal  Dielectric  Guide  Parallel  to  the  z- 
Axta.  The  rroaa-aertional  area  is 
denoted  by  A.  ita  bounding  curve  by  L, 
the  outward  normal  to  the  cylinder  ta 
n,  and  the  unit  vector  tangent  to  1.  la  t. 
The  poatttve  i-axta  la  out  of  the  plane" 
of  the  paper 


For  aolutiona  that  rorreapond  to  wave  propagation  along  the  guide,  it  will  be 
aaaumed  that  the  only  /-dependence  iae,kg7. 


(1) 


where  £  ta  a  vector  in  the  plane  perpendicular  to  the  t-axia  and  k^  ia  the  propaga¬ 
tion  conatant  in  the  1  -direction.  It  followa  that  E<£>  and  B<£)  aatiafy 


n 


t 


k 


o 


<fi> 


2  2 

Similar  formulae  hold  outside  the  guide,  with  t  j  and  >  ^  replaced  by  and  — ij. 
respectively. 

We  want  to  derive  tntegral  equation*  for  E?  and  It  and  want  thene  equation* 
to  incorporate  the  boundary  conditions,  which  are  that  the  tangential  components 
of  both  F  and  It  be  continuous  scross  the  boundary  curve  L.  These  tangential 
component*  are  the  longitudinal  Ey  and  and  the  transverse  and  B(,  To  get 
the  required  integral  representation  we  follow  a  procedure  very  similar  to  that  in 
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a  previous  paper  by  Eygea.  h  It  is  easy  to  find  an  integral  representation  for  E^ 
which  satisfies  Eqs.  (2)  and  (3)  and  also  guarantees  tliat  be  cor'inuous.  Thta 

is 


Et<P> 


*,  <£.£’>  dA' 


(7) 


where  the  two-dimensional  CJreen  funrtlon  satisfies 


IV*  -  ll)  g  (p.£'» 
i  *  ^2 


6(p  -  p’l  . 


(10 


There  is,  of  course,  a  similar  equation,  with  similar  propert‘es.  for  B  .  How¬ 
ever.  these  equations  do  not  satisfy  the  conditions  that  E(  and  B(  be  continuous. 

The  reason  they  do  not  is  that  the  condition  of  the  continuity  of  E(  and  B(  can  be 
rephrased  as  a  condition  on  the  discontinuity  across  the  boundary  of  3E^ /dn  and 

3n  (jump  condition),  whereas  Eq.  (7)  implies  that  E^  is  continuous,  and  the 
similar  equation  implies  that  B f  is  continuous.  But  Eq.  (7)  can  be  modified,  much 
as  in  the  previous  work  on  scattering,  bv  the  addition  to  It  of  certain  line  integrals 
that  guarantee  that  the  jump  condition  be  satisfied,  and  we  shall  do  just  this. 

To  get  the  jump  conditions  we  use  Eq.  (5)  for  E^  just  inside  the  boundary 
(subscript  -)  and  equate  it  to  its  counterpart  for  E^  Just  outside  the  boundary 
(subscript  •)  to  find 


Now  since  E^  is  continuous  across  the  boundary  so  is  3Ey  3t  so  we  do  not  need 

subscripts  •  and  -  on  it.  Thus.  Eq.  (9)  becomes  a  condition  on  the  discontinuity 

of  the  normal  derivative  of  B  ; 

z 


(10) 


In  much  the  same  way,  the  continuity  condition  on  B{  becomes  one  on  the  discon¬ 
tinuity  of  the  normal  derivative  of  E^: 


3 
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To  exploit  theae  reaulta,  we  recall  aorae  properiiea  of  the  following  line 
integral  over  the  boundary,  a  function  of  £  we  dealgnate  by  f(p): 


/s'* 


£')  <>(£•)  dl.' 


where  o  la  for  the  moment  an  arbitrary  function.  It  la  well  known  that  t(£)  la 
conttnuoua  acroea  the  boundary  but  haa  diacont tnuoua  normal  derivativea 


it  .  it 

in  ,  dn 


where  ^  on  the  right  hand  aide  of  Eq.  (13)  la.  of  courae,  on  the  boundary.  Now 
we  add  to  the  provtaiotial  Eq.  (7)  for  Kf  a  line  integral  like  Kq.  (12)  but  with  o(p) 
taken  to  be  the  function  on  the  right  hand  aide  of  Eq.  (11).  to  get: 


, p ')  dl.'  . 


Ef(*)  -  <>J  •  >j)  t  <£.£*»  dA‘  -  j-4  J  /(Cjtj  ‘  «2  ‘*1)  ^ 

A  *  2 ’I  1  l 

k  ,  ,  \ 

»  -jjl  (>j  ♦  ij)  -5 jt*  >  g^  (p.p*)  dl.'  (14) 

o  J  2"  " 

Thta  ta  now  the  final  deaired  equation  for  E  .  It  atill  aatlaftea  Eqa.  (2)  and  (3), 
guarantee#  the  continuity  of  E^,  and.  In  addition,  by  virtue  of  the  line  integral  in 
it  and  Eq.  (13).  now  aattafiea  the  jump  condition  Eq.  (11),  l.e.,  the  condition  on 
the  continuity  of  1^.  The  equation  for  Bf.  derived  alnillarly,  1a 

<£.£'*  dA’  — r 1  /  (SH  5^) 

A  J  ”*l  I.  '  °  ' 

g>  <£.g/>dl.'  .  (15) 

For  many  purpoaea  It  1a  convenient  to  tranaform  Eqa.  (14)  and  (IS)  by  replac* 
tng  the  area  integrala  with  line  integrate.  Thua,  In  Eq.  (14)  In  the  Integral  over  A 
we  can  uae  *£,')*  ^  g^  and  •  Ef  and  rewrite  It  aa 


we  can  uae 

followa: 
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*  «>?  *  >5'  /  E,  dA‘  *  /,Me  -  e'>  *  *2  c,  >  E*  ♦  /' ^  E« 

»  *  a  *  A  * 


*  E  » 

> 

.here  we  have  uaed  Green’a  eacmtd  identity  in  deriving  the  laat  equation  from  tta 
predecewaor.  With  Kq.  (16)  in  Kq.  (141  *e  find  after  cancelling  *<«me  it-nr  a 


/ 


<*K 


S,  liT 


-  E_ 


dn 


di. 


(p  tnatde  guide) 


(16) 


(  4,'tj  (  «.|  ,  “g  2 

0  ‘  J  Tr1  dl  *  '  iH  J  %  ^  d,  "77  K  (>»  * 

L  21!  1.  *  2M  ° 


/ 


w 


dl 


(17) 


In  much  the  lame  way.  Eq.  (IS)  hecomea 


0  •  J  Hr  l^2  <1‘-'  -“2  J  *n.  *7,  dl 

l  '  1  1  * 


1  1 


k  l  2 

O  V  S  J 


/ 

t. 


4E 

«y2  d*  ' 


(IS) 


Theae  are  the  haatc  integral  repi  rarntatione  of  the  theory. 


1  MJl 1%  MIST  SETS  t»»  U>»  \B  »QI  STIO>S 

In  thia  aectton  the  pair  of  integral  Kqa.  (171  and  (18).  ta  reduced  to  an  equiva¬ 
lent  pair  of  aeta  of  linear  algebraic  equal mna  b>  evaluating  them  at  amall  p,  much 
aa  waa  done  in  prrvioua  related  work  of  Evgra.  We  aaaume  that  p  la  much  leaa 
than  the  amalleat  value  of  p'  that  appeara  in  the  mtegranda  of  Kqa.  (17)  and  (18). 
W>  can  aaaure  that  thia  ta  poaaible  by  chnoalng  the  origin  of  the  coordinate  avatem 
in  Figure  1  appropriately.  Then,  for  uae  in  Kqa.  (17)  and  (18)  we  have  the  Green 
function  evpanaion 


IS 


Wf  ahall  uiumr  itvat  the  normal  n  ia  expreaaed  ut  cylindrical  coordtnatea, 

ri  *  p  n  ♦  i  n  (25) 

-  *■  *>  * 

although,  it  can  be  advantageoua,  particularly  with  a  cylinder  whoae  croas  aectlon 
ta  a  polygon,  to  expreaa  n  in  ('artealan  coordtnatea.  In  cylindrical  coordlnatea, 
for  a  contour  1.  that  ia  defined  by  a  pi*).  one  )iaa: 


'  '  a  - 

Theae  reaulta  are  put  into  Kqa.  (211  and  (22)  and  the  integration  ia  converted  to 
one  over  d#’  via  dl  '  (d!.'  da'lda',  whereupon  (hear  equal  Iona  become: 
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I 


-j  ]  7  ru'-n*  11,(1. 2p'l  VI/*>,P,>j|r 

V°  L 


The  homogeneous  set  of  Kq*.  (28)  and  (29)  ha*  solution*  only  if  the  determinant 
vanishes.  then  the  various  root*  of  this  vanishing  determinant  yield  the  propaga¬ 
tion  constants  of  the  different  mode*  of  the  guide. 

Eq*.  (28)  and  (2!*)  can  he  applied  to  circular  guide*  (radius  a)  to  obtain  re  - 
suit*  that  are  of  course  not  new.  but  which  do  serve  an  a  useful  check  on  our  pro¬ 
cedure.  f  or  a  circular  guide  <V>  .*♦  0  and  all  the  off-diagonal  matrix  element* 

are  rero,  that  i». 


M(«  a)Hj  (1>2*)*#>  (31) 

N#»  -*  •*'*2  JJ  <•>,•>  Hf  <1  >2  a>*f||  (32) 

•  2*<k  /koMI  *  I*  >*»)!  J,  <-» ,  •)  ff|  •  (33) 

If  we  let  x  yj*  *n<*  V  *  "*  j  »  tn  Eq*.  (31)  to  (33),  substitute  these  expressions 
into  Eqs.  (28)  and  (2*»).  and  then  eliminate  the  A(  and  I).  between  these  two  equa¬ 
tions,  we  get  the  standard  dispersion  relation  (Marcuse*); 


18 


(34) 


l(.  jk*  y)(-j;«y)  Jf(y))  -  («  2  k2/*MI}(t  x)/H#U  x))l  * 

|JJ(y)  y  j#(y)  -  IIJU  *)/«»!,  (lx))  I2  k2(y2  -  x2)2  /x4  y4  . 

*.  nm.nu  H»rK»»NT\noNs  for  weakly  (.mint.  hh).k.v  ■,  •  >2 


An  Important  practical  case  occurs  when  nj  •  so  that  the  fiber  is  weakly 

guiding.'  )'or  circular  flbera  Glogr'  lisa  shown  how  the  conventional  treatment 

can  be  simplified.  Mia  essential  point  is  that  the  longitudinal  components  Ef  and 

II  become  small  compared  to  the  transverse  components  K  ,  B  and  in  fact 
t  •  "  \  1  1 

1  ):  •  a  n.  -  n...  the  same  for  B  B  .  A  similar  result  holds  for  the  modes  of 

t  i  1  \  i  t  4 

the  slab  guide.  It  is  .  tear,  and  we  shall  verify  this  below,  that  Gloge'S  arguments 

cun  be  also  applied  to  weakly  guiding  fibers  of  arbitrary  cross -sectional  shape. 

It  then  becomes  more  convenient  to  abandon  the  treatment  of  Section  3  in  terms  of 

E  and  B^  and  to  deal  directly  with  the  nonvanishing  transverse  Helds  E^  and  B^. 

We  want  then  to  obtain  integral  representations  for  E  and  B  analogous  to 

those  of  Section  3  for  )  and  B  .  To  set  up  these  latter  integral  representations. 

«e  had  to  take  into  account  the  continuity  properties  of  their  normal  derivatives 

across  the  boundary.  Similarly,  in  setting  up  the  integral  representations  for  E 

and  B  *e  must  again  Inquire  into  the  possibility  of  boundary  discontinuities.  To 

begin,  we  observe  that  the  tangential  components  of  E  and  B  must  lie  continuous 

serosa  the  boundsry.  What  about  the  normal  derivatives  of  these  qusntittes  ’  To 

studv  this  we  set  up  an  orthogonal  n.  t  coordinate  system  with  Its  origin  st  some 

point  on  the  boundary,  wherein  the  n-axls  is  slong  the  direction  of  the  normal  n, 

at  that  point,  and  the  t-axls  i«  along  the  direction  of  the  tangent  vector  Then, 

in  terms  of  these  variables  the  t  -component  of  C  >  B  just  outside  the  guide  is 

<r  V  B)  (AB  An  -  AB  At)  IkB  -  0.  We  can  write  a  similar  equation  just 
—  —  *  t  n  ♦  i 

inside  the  boundary  and  hence  conclude 


(AB.  AB  \  /AB.  AB  \ 

- 1-  - &  [ - 1.  - D  ]  (3  5) 

An  At  I  l  An  At  I 

Now  r  -  B  0  which  means  that  B  is  continuous  across  the  boundary,  and,  hence. 
—  —  n 

that  its  tangential  derivative  AB  At  1"  also  continuous.  Then  Eq.  (35)  becomes 

n 


(36) 


IP 


In  short,  the  normal  derivative  of  the  tangential  component  ia  continuoua  acroaa 
the  boundary.  Consider  the  rectangular  component  a  B  and  B  .  If  B  and  B  and 
their  normal  derlvativea  are  continuoua  acroaa  the  boundary,  then  the  tangential 
component  of  B  (a  linear  combination  of  B^  and  IO  will  alao  have  theae  con¬ 
tinuity  propertiea.  Arguments  aimtlar  to  the  ab<ivr  alao  apply  to  E^  and  Ey. 


Let  •  be  either  of  tl*.  Bx<  By. 


C'onalder  Eq.  (7)  vnth  Ef  replaced  by  ♦: 


•  <£>  ■  -  <1, 


7?* 


f 


(£,£*)  dA 1 


(37) 


As  we  have  seen,  this  guarantees  that  Eqs.  (2)  and  (3)  are  satisfied  and  that  t 
and  d#  dn  will  be  continuous  across  the  boundary  (Eyges  and  Nelson*®).  If  thia 
equation  is  aattafled,  we  can  take  •  to  be  either  B^  or  B  .  or  n  linear  combination 
of  the  two  of  them.  This  linearly  polarlred  mode  *111  then  propagate  down  the 
guide  of  arbitrary  crosa  section  *tth  propagation  constants  and  modal  character¬ 
istics  that  are  governed  by  Eq.  (37)  and  with  an  E  field  that  is  perpendicular  to  the 
B  field,  as  determined  by  Maxwell’s  equations.  Note  that  Eq.  (37)  can  be  converted 
to  a  useful  line -integral  form  by  the  steps  that  led  from  Kq.  <lf>)  to  Eq.  (17).  It  be¬ 
comes,  for  p  inside  the  guide. 


/ 


-  dn’ 


-it 

i*n’ 


dl.’ 


(38) 


It  is  worth  observing  that  Eq.  (38)  can  be  derived  in  an  alternate  way  that  ia  re¬ 
lated  to  the  so-called  Extended  Boundary  Condition  method  (EBC)  which  has  re¬ 
cently  been  exploited  for  scattering  problems  (Waterman**).  Basic  to  this  method 

is  an  integral  representation  (tluygen’a  principle)  for  an  amplitude  #  which  satts- 
2  2 

flea  (V  -  0  in  the  exterior  region,  namely 


•  (£) 


f  L  .  a*  \ 
j  y*  **>•  *  *>2/ 


dl.’ 


£  outside 


(3P) 


It  is  well  known  that  if  this  equation  la  evaluated  in  the  interior  (canonical  proce¬ 
dure  In  the  EBC  method).  It  yields  rero  on  the  left  hand  side: 

10.  Eyges.  I.,  and  Nelson,  A.E.  (1976)  Ann.  Phys.  100:37. 

11.  Waterman.  P.C.  (1969)  J.  Acoust.  Soc.  Am.  45:1417  and  (1971)  Phys.  Rev. 

D3:B25.  - 


Hut  for  thr  present  problem  we  have  _  and  i)€  9n|  4  •  84  8n|  _  and  when 

thear  taro  equation*  are  aubatituted  into  the  above,  the  baaic  Kq.  (38)  la  reproduced. 


1  xlMMITKl  COMSIDf  RATIONS 

Kither  of  the  integr  al  repreaentationa  (37)  or  (38)  can  noar  b«  converted  to  a 
set  of  homogeneous  linear  equation*  and  a  corresponding  determinantal  equation 

I  HA 

b>  inserting  m  them  an  expansion  in  the  interior  solutions  J^lt  jp)  c  .  In  fact, 
it  is  more  convenient  to  uae  trigonometric  functions,  so  «e  shall  assume  aa  the 
general  interior  expansion  for  4  for  a  guide  of  arbitrary  rroae-aectional  shape: 

* 

•  (£)  y]  J.<7tp)  (C.  cos  (*4>  *  PM  sin  (ad)  I  (40) 

a  -o 

h  q.  (40)  is  in  fai  t  too  general  in  that  almost  all  practical  guides  have  symmetry  of 
one  kind  or  another.  It  is  then  useful,  indeed  almost  imperative,  to  sort  out  In 
advance  the  limitations  that  such  symmetry  implies  for  the  coefficients  in  Kq.  (40). 

In  this  paper,  for  example,  we  shall  analyse,  among  others,  guides  of  rec¬ 
tangular  or  elliptical  cross  section.  If  the  axes  of  such  an  ellipse  or  rectangle 
coincide  with  the  x-  and  v-axes,  then  by  usual  parity  arguments  of  quantum 
mechanics  the  amplitude  4  must  be  either  symmetric  or  antisymmetric  about  the 
x-axis.  This  means  tliat  4  can  be  expanded  in  terms  of  eit her  sines  or  cosines. 
Similarly,  reflection  symmetry  about  the  y-axis  ran  be  used  to  shot*  that  both  the 
sine  and  the  .  sine  expansions  split  Into  two  other  rases,  in  one  of  which  s  runs 
over  even  values  and  In  the  other  over  odd  values.  In  short,  there  are  for  ellip¬ 
tical  or  rectangular  cross  sections  (or  any  guide  with  similar  twofold  reflection 
symmetry)  four  t\}>e*  of  expansions,  labelled  by  R  (for  rectangle)  with  superscripts 
I  through  IV: 


L 


C’s  jp)  cos  <s4) 


(41s) 
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C«  Ja(,>  lp)  co“  <M> 


(41b) 


H  :  t(£) 


*  1.  3,  5. . . 


H,U:  •(£) 


in  <a«) 


H  2.  4.  6. . , 


R,V:  #(p)  Ds  VV*  *ln  (,*#) 


a  I.  3.  5,  .  . . 

If.  in  .«  ddition  to  the  twofold  symmetry  in  the  jierpendlcular  axes  of  the  rec¬ 
tangle  or  ellipse,  the  guide  has  the  symmetry  of  the  square  (reflection  symmetry 
with  respect  to  diagonals),  the  classifications  of  Kq.  (41)  are  modified  In  a  way 
that  cannot  be  derived  by  quite  so  simple  reasoning  as  the  parity  arguments  above. 
The  results  are.  however,  a  straightforward  consequence  of  group  theory  and  of 
the  fact  that  the  symmetry  of  the  square  is  described  by  the  group  t '  which  has 
(  ■Ur  onr-dtmrnai<  nal  trredu<  ible  representations  and  one  two-dimensional  one. 
This  leads  to  the  conclusion  that  there  are  four  expansions,  namely,  Kqs.  <42a-d), 
which  lead  to  four  different  (nondegenerate)  sets  of  propagation  constants,  and 
two  expansions  (Kqs.  <42e  and  0)  that  yield  that  same  set  of  propagation  constants 
(degenerate  modes).  These  are 


•  <£>  ^  t  ^  J^ly  jp)  cos 


*  0.  4,  8. ,  . 

♦<£>  2 

s  2.  «.  10.  . 


C's  J#(>  .P >  cos  (s$> 


•<e»  2 


**r* 


s  2.  0.  10.  . . 


,p-»*  £  °.V*l 


p)  sin  (*«) 


s  4.  8.  12... 
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c*  Js<,> lp*  co* 


<42e> 


H"  ♦<e)  £ 

s*  1.  3,  5. . . 


HIV  •(£) 


E 


Dg  (•#)  . 


(42f> 


However,  in  order  to  relate  the  mode*  of  the  square  to  thetr  counterparts  of  the 
rectangle,  we  have  kept  the  nu>de  notation  introduced  in  Eq.  (411.  Thus,  in 
Eq.  (42)  R*  at  til  refer*  to  "cos,  even  s”  modes  and  It'*,  to  "cos,  odd  s"  modes, 
etc. 


6.  r.yi  4T1«ns  eo*  the  weakly  c;i  iihm.  case 

Given  the  expansions  of  Eq.  (41)  or  (42),  we  can  insert  any  of  them  into  the 
integral  representations  (37)  or  (38)  and  use  the  expansion  of  the  Green  function 
to  obtain  a  set  of  linear  equations  and  hence  a  determinants!  equation  for  the  modal 
propagation  constants. 

One  question  of  strategy  is  now  whether  to  chooae  the  line  integral  represen¬ 
tation  (38)  or  area  integral  representation  (37).  The  former  choice  might  appear 
to  he  the  obvious  one  since  the  line  integral  matrix  elements  it  generates  are  one* 
dimensional,  and  hence  ostensibly  simpler  than  the  two-dimensional  integrals  that 
Eq.  (37)  would  yield.  In  fact,  they  are  generally  simpler  and  we  shall  choose  to 
work  with  Eq.  (38),  but  Eq.  (37)  should  not  be  set  aside  completely;  It  Is  particu¬ 
larly  useful  for  guides  that  are  not  too  far  from  circular.  In  this  esse  it  lends 

12 

Itself  to  a  convenient  perturbat Ion  analysis  (Evges  )  in  which  the  integration  la 
broken  up  into  one  over  a  circle  (for  example,  inscribed  within  the  cross  section, 
or  of  the  same  area)  which  generates  diagonal  matrix  elements,  and  an  Integral 
over  the  remainder  area  which  yields  off-diagonal,  but  small,  matrix  elements. 

For  the  present  work  with  Eq.  (38)  we  shall  need  the  expansion  of  the  Green 
function  in  sines  and  cosines: 


«■>  't-t''  *  -J  ^  *,  H|(iy2p')J#  <y2p)|coS  14  cos  14'  *  sin  14  sin  / ♦*  1  . 


I  *o 

p  <  p' 


(43) 


12.  Evges,  L.  (1978)  Appl.  Opt.  L7:1673. 
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where  «  ^  the  Neumann  factor,  is  1  if  I  0  ami  2  otherw  ise.  The  details  involved 
in  substituting  Kq.  (41)  or  (42)  into  Kq.  < H H )  are  similar  to  those  spelled  out  in 
Section  1.  The  cosine  expansions  ultimately  yield  the  system  of  linear  homogene¬ 
ous  algebraic  equations: 

^  ' C  G.  0  (44) 

s  Is 

a 

whereas,  the  extensions  involving  sines  eventually  produce 

EOT,  0  .  (45) 

s  la 

a 


The  indices  I  and  s  in  the  above  rqunt  ions  are  either  N>th  even  or  both  odd  and 
take  on  the  values  spei  ifted  in  Tq.  (41)  or  (42).  G(  ^  and  T/n  which  are  line  inte¬ 
grals  taken  about  the  perimeter  of  the  guide,  are  given  by: 

2» 

r.lH  /  dp  |Kf#(p')  cos  Ip*  coa  sp'  »  l  (>(p\p  ) 

a  O 

(a  cos  (s'  Bln  a#'  -  I  Bin  1$'  cob  sp'}|  (46) 

2» 

Tfn  f  dp’|Ff  Jfi')  Bin  »♦•  ain  a#'  -  L(-(p',*') 


(a  aln  Ip'  cob  bb'  -  I  coa  14'  aln  ep')l  (47) 

with 

F(  #<p')  c|t>,p1H|(t>2P*)j;(>1p')  -  tjP’HJ  <1'.2p,)J(|(y1i»')|  (46) 

,  la<p’-#')  rl  H» l^/p'  •  (i9) 

The  related  determinants!  equations  are 

det  (G.  )  0  (50) 

la 

del  (T  )  -  0  .  (51) 


24 


One  >'in  easily  recover  the  familiar  dispersion  relation  for  a  weakly  guiding 
circular  guide  of  radiua  a  from  Kqs.  (44)  to  (49).  Cl  become#  diagonal  (and  iden¬ 
tical  to  T)  becauae  l  0  and  K(<  come#  outside  the  Integral.  Then  the  condition 

yil  1  «>»,  <1  >2  “>J,|  <■),»>*  (1,  a)  •  o  (52) 

hold#  for  all  I  and.  when  satisfied,  give#  the  modes  of  the  circular  guide. 


:  conn  r \no>  \i.  TM.HMyi  rs 

In  this  section  we  outline  some  of  the  detail#  of  our  computational  technique#, 
l  or  the  work  of  this  paper,  Kqa.  (48)  and  (49)  are  im>dtfied  In  three  way#.  First, 
the  lies  sc  l  function  derivative#  are  eliminated  by  u#e  of  the  general  recursion 
relation 

Cj(x)  <■.<*)- ~  Cf<*>  •  (53) 

Then  the  Hankel  function#  are  replaced  by  Modified  Bessel  functions  of  the  second 
kind,  using 

H#(!  x)  i  ,*<»*>>  ^(x)  .  (54) 

Finally,  the  arguments  i,  and  •>..  are  expre##ed  in  term#  of  a  normalised  props - 

2  1  *  2  2 

gallon  constant,  P  ,  and  a  normalized  binding  constant,  ;1.  P  ,  following  CJoell, 

has  tern  normalized  so  that  it  varies  between  0  and  1  only,  while  ;1  characterizes 
the  refractive  index  dlfferem  e,  the  smaller  guide  dimension  and  operating  fre¬ 
quency,  hence  the  extent  to  which  the  field  exist#  outside  the  guide.  These  quan¬ 
tities  are  defined  to  be 


I* 


2 


bk 


(nt 


(55) 


<5fi> 
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* here  b  m  the  aeminunor  *xi»  of  the  cross -sect tonal  shape.  *  One  should  note 
that  with  these  definitions  specific  values  of  b,  k  ,  nJt  and  n.,  are  not  required. 
<  me  then  obtains 


ijP  *3  ;  tjd  '  P|lp/b  . 


With  thes«  niiHlifU  allons,  Fqs.  <4H)  and  (49)  t  an  be  rew  ritten  as 


v,  (H*  «7TP)  vr 

,PJ_  (^.  VT7)  (^)J 

4.-»-"k,  (^)x. 


To  extract  numerical  results,  four  things  must  be  done,  First,  ptp)  must  be 
specified,  second,  the  matrix  elements  C.(  r  (or  T(  )  must  be  computed;  third,  the 
roots  of  the  determinants!  equation  must  be  found,  and  finally,  for  each  root,  the 
coefficients  (‘b  (or  l)^)  are  to  be  found  to  within  a  normalising  constant  so  that  the 
field  configuration  can  be  determined  from  Fq.  (41)  or  (42). 

To  calculate  the  line  integrals  that  comprise  the  matrix  elements  there  are 
two  computational  possibilities,  (me  ran  express  the  integrands  In  rectangular 
coordinates  and  integrate  over  dl.  directly,  or  one  can  convert  the  line  integrals 
to  angular  ones  by  writing  dl  (dl  d^l  d4  and  integrate  from  ;ero  to  2*.  The 
latter  procedure  is  clearly  desirable  for  curved  boundaries  (for  example,  ellipses) 
for  which  p  and  Sp  '<*♦  are  continuous  functions  of  4.  but  for  rectangular  cross 
sections  the  direct  ini  ^ration  seems  more  straightforward.  For  the  sake  of 
generality  in  the  computer  programming,  however,  we  have  calculated  the  matrix 
elements  even  for  rectangular  boundaries  by  angular  Integration,  using  the  function 


p(#)  •  b|(cos2  #  ^)N  *  (sin2  #)NI  ‘l  2N 


With  this  definition.  3  turns  out  to  be  one-half  of  Caoell’s  corresponding  quantity, 

B. 


— w— 


fur  the  perimeter.  This  uivm  a  figure  with  a  short  axis  (along  p  i»  21  of  length 
2b  and  a  long  axis  (along  «  0.  »)  equal  to  2Rb,  where  R  is  the  aspect  ratio.  In 

general.  for  arbitrary  H.  Eq.  ((>0)  describes  a  rectangle  if  the  |>arameter  N  is 
i  hoseri  to  be  infinity  (Allard*  *  and  Gardner**)  and  an  ellipse  If  N  1.  Sha pes 
with  l  <  N  <  «■  are  called  " suj>ei  ellipses  '  (Gardner**)  and  can  be  regarded  as  a 
continuous  deformation  of  an  ellipse  Into  a  rectangle  as  N  increases.  Shapes  with 
0.  5  *  N  *-  l  are  called  '  subellipses"  (Gardner  ’)  and  can  be  regarded  as  a  contin¬ 
uous  deformation  of  an  ellipse  into  a  parallelogram.  \  0.  5  gives  a  parallelogram, 

and  if  0  *  \  (  0.  S  one  obtains  figures  w  ith  four  cusps,  tending  toward  crossed 
straight  lines  as  N  -0.  Whenever  R  I  in  the  above  analysis,  the  rectangle  and 
the  ellipse  become  the  square  and  the  circle,  respectively. 

In  figure  2  we  show  the  first  quadrant  of  some  of  the  shapes  obtainable  for 
various  \  with  R  1.  Hegardless  of  the  values  of  N  or  R  in  Eq.  ( (10),  reflection 
symmetry  about  the  x-  and  v-axes  produt  es  a  closed  figure  throughout  the  four 
quadrants.  Thus,  use  of  Eq.  (*0)  assures  continuity  of  dp  d#  over  the  entire 
range  of  o.  thereby  eliminating  the  need  of  subdividing  $  into  several  intervals 
whose  boundaries  depend  on  R. 

The  fractional  error  >>nr  makes  in  taking  Eq.  <*0)  as  an  approximation  to  a 
ret  tangle  is.  for  a  given  value  of  N,  independent  of  R.  The  error  in  p<$)  i»  worst 
at  the  corners  but  becomes  small  very  rapidly  for  values  of  $  awa\  from  the 
corners.  We  used  V  10,  which  gives  a  fractional  error  of  0.012  at  the  corners, 
0.002  one  degree  nwav,  ami  0.00025  two  degrees  away,  f  or  N  10  the  values  at 
the  same  three  angular  positions  are  0.014 .  0.020 ,  and  0.011  respectively. 

Hecause  of  symmetry,  the  integrals  (41  »  and  <47)  need  to  be  done  only  over 
the  interval  (0.  »/2).  W  ith  the  transformation  t  -cm  they  become  suscept - 
ible  to  Chebyshev  integration  iMildebrand  )  which  takes  equally  spaced  points  in 
the  range  (-1,  1)  with  equal  weights.  Generally.  0  points  were  sufficient  to  obtain 
the  required  a>  curacy,  although  t  hecks  were  made  with  many  more. 

In  order  to  satisfy  Eq.  ( SO)  or  (31),  the  determinants  were  truncated  and  then 

2 

evaluated  at  equallv  spaced  values  of  J’  .  The  rinds  were  located  approximately 

from  the  table  of  these  results,  and  a  half-interval  search  (Hornbeck*®)  was  uaed 

2 

to  pin  down  the  values  of  1’  to  four  significant  figurea.  The  number  of  terms  to 
retain  was  determined  empirically,  with  two  terms  being  sufficient  in  many  cases. 
For  a  given  mode  of  a  given  guide,  the  number  of  terms  needed  for  a  certain 
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)  igure  2.  Plot  of  the  function  p (o >  (liven  in  Kquat  ion  (fiO)  for 
0  -  e  -  »  2  Only.  A  J,  and  Various  Values  of  N 

2 

absolute  irrur»'  V  in  I’  dofu  not  depend  on  Ilmevor,  in  general,  different 

require  different  numhrrt  of  term*  to  jive  rrmilli  with  the  nmf  accuracy. 
The  computation*  acre  carried  out  on  a  C'tX'-fiSOO  computer.  The  time  required 
la  nearlv  proportional  to  the  square  of  the  dimension  of  the  determinant,  for  a 
g’ven  numler  of  roots  and  accuracy.  Thus,  a  three-term  truncation  procedure 
might  use  .10  sec  to  loi  ate  three  roots  to  within  0,0001.  t'aing  four  or  more  terms, 
it  is  usually  possible  to  search  a  limited  region  near  a  root  and  thus  save  time. 

a  RK51  ITS  AM)  COW  ARISON  WITH  OTUKR  DORK 

We  have  investigated  the  properties  of  weakly  guiding  waveguides  having  all  of 
the  various  cross  sections  described  in  Section  7  and  aspect  ratios  ranging  between 
I  and  10.  In  this  section  we  present  some  of  our  reaulta,  to  demonstrate  the  util¬ 
ity  of  the  method  and  make  comparisons  with  previously  published  work  for  the 
more  common  shapes. 


d  I  IrrtMfllUi  Otws  SkImh 

For  the  rectangular  wavr  guide*  having  aspect  ratio*  t  -  ft  t  2  we  have 

located  more  than  two  doxen  mode*  in  the  range  0  «.  fl  *  4.  he  curve*  for  the 

2 

propagatu.y  constant*  (P  ver*u*  ,i)  corresponding  to  approximately  the  flrat 
dozen  mode*  for  each  aspect  ratio  are  plotted  in  Flgurea  3,  4.  5,  and  6  for  if 
value*  up  to  2.  5.  Fach  of  these  curve*  correspond*  to  one  of  four  type*  of  field 

function  when  ft  >  1  (see  Eq.  (411)  or  one  of  six  type*  when  ft  1  (wee  Eq.  (42)). 

I  IV 

The*e  field  function*  are  designated  H  through  It  The  mode  cor  responding  to 

•he  n-th  root  of  the  X-th  type  of  determinant  (see  Fq.  (SO)  or  (SI))  i*  designated 

X  1 

H'  (Thu*,  the  first  (noitcutoff)  mode  i*  labelled  H,  because  It  involves  a  cosine 
n  i 

expansion  with  even  indices.  )  The  ordering  within  a  family  is  governed  bv  the 
value  of  fl  at  which  cutoff  occur*,  n  I  designat  ing  the  nu«dc  which  cut*  off  at  the 

smallest  .t.  Since  for  ft  »  1  certain  mode*  cross  over”  (for  example,  see  M*  and 

I  ■ 

in  Figure  ti).  this  choice  must  be  kept  in  mind  when  considering  modes  at 

value*  of  .i  above  the  crossover  point  (for  example,  see  the  values  of  ;i  tor  these 

same  two  modes  In  Table  4). 

As  was  noted  in  Section  5,  the  special  symmetry  of  the  square  cause*  the  even 

expansion*  to  split  into  two  series  each.  Thu*,  consider  the  twro  cosine  expansions 

of  Eq*.  (42*  and  b):  The  first  rm>de  determined  using  Eq.  (42a)  is  the  first  mem* 

ber  of  ",e  M!  family  (that  is.  It')  and  its  field  expansion  is  given  as  C  J  4  C .J  . 

1  o  o  4  4 

cos  4«  4 .  Whereas,  the  first  mode  determined  using  Eq.  (42b)  turns  out  to 

be  the  third  member  of  the  same  family  (hence,  Hj)  whose  field  exfuinslon  has  the 

form  t  COS  2*  •  cos^s  4 .  However,  for  the  rectangular  cross 

section*  (ft  ■>  l),  the  field  expansions  for  all  even  cosine  modes  take  the  form 

C  J  4  C_J  -  cos  2*  4  C‘ .J  .  cos  44  4 . 

o  o  l  l  4  * 

A  crucial  question  underlying  our  approarh  ia:  How  many  terms  must  be  kept 
in  the  determinant*!  Fq*.  (50)  and  (51)  in  order  for  the  root*  to  converge  to  de¬ 
finite  values'*  U'e  now  address  this  question. 

The  convergence  of  the  fl!  expansion  la  rapid  for  both  the  square  and  rectan- 

1  2 

gular  guides.  For  the  square  guide.  Table  1  give*  the  values  of  P  obtained  for 
this  mode  after  truncation  of  F.q.  (SO)  at  1.  2.  or  3  terms  for  selected  values  of  fl. 
Clearly,  the  tso-tertn  truncation  I*  only  a  small  improvement.  In  terms  of  locating 
the  roots,  on  the  single -term  result,  which  is  obtained  by  setting  C.qq  *  0  In 
Eq.  (50).  Moreover,  the  extra  work  involved  in  taking  three  terms  leada  to  no 
gain  whatsoever.  In  view  of  this,  it  is  reassuring  that  the  ratios  of  the  second 
term  of  the  field  expansion  to  the  first  (C'^J^ ,  C'qJq).  and  the  third  to  the  first 
(CgJg  CpJg)  Show  that  the  term  is  dominant  for  this  mode.  The  term  ratios 
listed  In  the  5th  and  fith  column*  of  Table  1  were  evaluated  on  the  perimeter  at  the 
middle  of  one  side  (th*(  is.  p  b.  4  0).  Since  and  become  smaller  as  p 
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Figure  4.  Mode*  of  the  Weakly  Guiding  Rectangular  Guide  with  ft 


1.0 


p* 


Q« 


0* 


o.« 


a? 


o 


050 


•00 


0 


I  50 


7  00  J 50 


l  igurr  s.  \|n|r«  of  thr  Weakly  Guiding  Rectangular  Guide  with  R  1.5 


f  igure  *i.  Modes  of  the  Weakly  Guiding  Rectangular  Guide  with  R  2.  Square 
point*  •  ere  taken  from  Figure  17  of  Goell.  *  X’*  repreaent  result*  of  a  per¬ 
turbation  calculation  for  the  Rl  mode 
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Table  1.  Hrsulti  for  the  Dominant  (Noncutuff)  K,  Mode  of  the  Weakly  Guiding 
Square  Guide.  Luted  value*  of  P-  were  obtained  after  a  1-,  2-,  or  3-term 
truncation  of  the  field  expansion.  Term  ratio*  derived  from  the  3-term  aolutlon 
give  G4J4  co»4*  C'oJq  and  cm  84  C'oJo  11  P  b.  $  0.  The  field  intenaity 

at  p  t  0  la  al*o  given.  The  intensity  at  p  0  i*  1 


P2 

Term  Hattos  at 
Perimeter  ($  01 

Intensity  at 
Perimeter 
<0  0) 

3 

1  Term 

2  Terms 

1  Terms 

‘  4J4  ‘Vo 

‘V«  ‘Vo 

m 

0.  9511 

0.  9546 

0.  9546 

-0.  314 

0.  0006 

0.  023 

m 

0.  9329 

0.  9371 

0. 9372 

-0.  266 

0. 0009 

0.032 

2.0 

0.  9023 

0.  9072 

0.  9072 

-0. 213 

0. 000  :> 

0.046 

1.  5 

0.  8442 

0.  8497 

0. 8497 

a 

O 

V* 

UJ 

0.  0005 

0.  076 

1.0 

0. 7122 

0.  7175 

0.  7175 

-0. 089 

0.  0004 

0.  143 

0.  5 

0.  3263 

0.  3290 

0.  3291 

-0. 029 

0.  0003 

0.  371 

decreases,  the  term  ratio*  are  generallc  greatest  at  the  perimeter  and  as  such 
constitute  the  worst  case  for  convergence.  (For  ,i  3  at  (p,*t  (0.  5b,  01  the 

term  ratio*  are  -0.007  and  *  10  ’,  respectively.  The  intensity  is  0.576.) 

Table  1  als.  list*  the  field  intensity  at  p  b,  c  0,  normalised  to  a  value  of 
1  at  the  center  of  the  guide.  For  the  more  tightly  bound  ra*cs,  (higher  ,1),  the 
field  drop*  off  to  n  few  per <  ent  a!  the  edge,  whereas  for  weakly  bound  case* 
(smaller  .0  the  intensity  is  still  appre<  sable  there,  indicating  that  a  greater  frac¬ 
tion  of  the  beam  1*  extending  beyond  the  guiding  core.  One  can  al*o  see  that  for 
those  cases  in  which  the  second  term  (C^.l^  C 'p.T ^1  makes  the  largest  [<ercent>  con¬ 
tribution  to  the  field,  the  total  field  i*  small,  whereas  for  weaklv  bound  <  ase*  the 
second  term  is  a  much  smaller  fraction  of  a  considerable  larger  field.  (Similarly, 
in  the  interior,  the  percent  contribution  of  the  sc.;  ond  term  drops  quickly  as  the 
field  Intensity  increases.  >  Viewing  the  guide  as  a  whole,  therefore,  all  terms 
bevond  the  first  one  mske  unimport  ant  rontribut  ions  to  the  total  intensity  for  the 
dominant  mode  of  the  square  guide.  The  same  holds  true  for  rectangular  cross 
sections. 

In  Figure  7  the  field  intensity  is  plotted  as  a  function  of  p  b  for  this  mode  both 
along  a  line  th;  .  ugh  the  renter  of  the  perimeter  (#  0°)  and  along  a  diagonal 

(«  45°). 

We  also  obtain  gi*od  convergence  of  the  field  extensions  for  the  remaining 
(rutofn  modes.  Table  2  compares  the  l -,  2-,  3-,  and  4-term  truncation  results 
of  F.q.  (50)  or  (51)  for  the  first  severs!  models  of  the  square.  Comparing  the 


FIELD  INTENSITY 


P/b 

Figure  7.  field  Intrn/tit \  as  a  Kuncl  ion  of  p  h  for  the  III  Mode  of  a  Square  (iuide 
with  ,1  3.  The  l*o  value*  of  #  rhuien  represent  the  extreme*  of  the  intensity 

distribution 


2-  and  4-term  truncations,  we  see  that  ont>  two  terms  are  required  for  an  accu- 
2 

racy  of  0.01  in  P  for  7  of  the  first  13  modes  of  the  square.  As  for  the  remaining 

modes,  three  term*  are  sufficient  to  bring  the  arruraev  to  within  0.01.  except  for 

the  degenerate  M - 17 1,'  mode.  which  is  accurate  to  within  0.02. 

Table  2.  which  also  list*  the  mode  designations,  specifies  in  it*  last  column 

which  term  is  dominant  in  the  field  expansion.  This  is  indicated  by  the  underlined 

value.  Thus,  we  see  that  for  most  of  the  modes  the  first  term  does  dominate,  and 

It  is  this  term  used  alone  in  Fq.  <30)  or  t  M  >  whirh  produres  the  1-term  truncation 

result  in  Table  2.  For  example.  (i__  0  was  used  for  the  R*  mode,  T~»  0  for 

III  III  ‘‘  ■’  **  II  iv 

'he  R"  .  T^  0  for  the  R"1  and  either  ('.jj  or  Tjj  0  for  the  degenerate  R" -Hj 

mode.  The  observation  stated  above  concerning  the  R *  mode  is  generally  applic¬ 
able  to  these  other  modes.  That  is,  for  mini  inodes  the  higher  terms  in  the 
expansion  are  comparable  in  magnitude  to  the  first  term  only  where  the  Reid 
Intensity  is  small;  otherwise,  thev  are  negligible. 


3  3 


2 

Table  2.  Values  of  the  Propagation  Constant,  P  ,  for  3  -  2  for  the  Ftrat  13  Modes 
of  a  Square  Guide,  Obtained  with  Determinants  up  to  Order  4.  Modes  are  ordered 
according  to  the  value  of  3  at  cutoff.  Our  mode  designation  and  Goell's  are  listed, 
and  the  dominant  term  in  the  field  expansion  is  indicated  by  an  underlined  expan* 
slon  index.  With  1-term  truncation,  the  dominant  term  only  is  used 


Mode 

Designs!  ion 

„2  . 

This 

3  at 
Cutoff 

r  ax  p  •  i 

Kx  pans  ion 
Indices 

Work 

Goell* 

1  Term 

2  Terms 

3  Terms 

4  Terms 

*5 

F«-V 

0 

0.  9023 

0.  9072 

0. 9072 

0. 9072 

0,  4.  8, . . . 

»!' 

0.  70 

0.  7  S4« 

0.  7614 

0.  7688 

0.  7 69 S 

1,  3.  5.  . . . 

«r 

K*.  y  / 

rl  2  J 

\j  3,  5, . . . 

H1” 

kx*  y 
r‘2  2 

1. 04 

0.  6107 

0.  6321 

0.  6323 

0.  6323 

2,  6,  10,  . . . 

R* 

K«.y 
r  3  1 

l.  10 

0.  5037 

0. S410 

0.  S4S7 

0.  S4  SB 

0,  4.  8 _ 

»; 

rx.  V 
f  1  3 

1.  17 

0.  S3 90 

0.  S446 

0.  S446 

0.  S446 

2,6.  10.  ... 

H*V 

1.41 

0.  3  383 

0.  37S4 

0.  3 8 SO 

0.4071 

1.  3,  5 _ 

1.3.  5 _ 

«? 

\ 

1. 62 

0.  2300 

0. 2230 

0.  2461 

0.  2447 

l  3.  5.  . . . 

"!v 

/ 

1,  3.  5 _ 

l.  74 

0.0784 

0.  IS08 

0.  1892 

0. 1902 

0.  4,  8.  . . . 

R»‘ 

1.81 

0.  0S9 1 

0. 1118 

0.  1240 

0. 1243 

2,  6.  10.  ..  . 

R»" 

1.  87 

0.  0914 

0.  1091 

0.  1089 

0.  1089 

4,  8.  12.  ... 

However,  there  are  some  modes  for  which  the  second  term  In  the  field  expan¬ 
sion  is  the  moat  important  one.  For  example,  the  field  Intensity  of  the  R^  mode 
of  the  square  guide  is  plotted  in  Figure  H.  It  shows  peaks  at  two  different  radial 
distances  from  the  center.  In  the  peak  near  the  origin,  the  CgJ^  term  turns  out  to 
he  the  dominant  one,  but  in  the  second  peak  the  first  term  ratio  has  values  ranging 
up  to  20  and  is  consistently  greater  than  1,  showing  that  the  C^J^  cos  44  term  is 
the  most  important  one  there.  For  this  mode  setting  only  G^  «  0  produces  the 
l -term  truncation  result  shown  in  Table  2.  Another  example  In  which  the  second 
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Figure  8.  Field  Intensity  an  a  Function  of  p  b  for  the  R^  Mode  of  the  Square 
Guide  with  J  .1  and  Various  Values  of  p.  When  p  >  0.4b  the  second  term 
(C^J^  roi  4#>  dominate*  the  field  expansion.  (See  text) 

IV  11 

term  predominate*  la  the  degenerate  R^  -R^  mode.  Here,  the  1 -term  truncation 
result  wa*  obtained  by  setting  either  Gjj  or  0.  In  such  cases  as  these.  It 

is  obvious  from  a  perusal  of  the  Table  2  result*  that  a  I -term  truncation  la  Inade¬ 
quate  for  proper  convergence  and  that  at  least  three  terms  In  the  expansion  are 
necessary. 

The  mode  R^  1  t*  an  example  tn  which  even  though  the  first  term  la  dominant, 

the  other  terms  are  still  making  a  significant  contribution.  From  Table  2  one  can 

see  that  here,  also,  three  terms  are  required  In  the  expansion. 

Tables  3  and  4  give  the  same  type  of  Information  for  the  R  •  2  rectangle  as 

Tables  1  and  2  did  for  the  square.  One  would  expect  the  convergence  at  a  given 

order  of  truncation  to  become  tncreaatngly  worse  for  higher  and  higher  modee. 

2 

This  I#  borne  out  in  Table  4.  Here  are  shown  the  effects  on  P  of  truncating  up 
to  seventh  order  for  the  first  13  modes  of  the  R  •  2  rectangular  guide  at  a  fixed 
value  of  B.  These  results  demonstrate  that  either  a  2-  or  3 -term  truncation  is 
sufficient  for  an  accuracy  of  0.01  or  less  in  P  for  7  of  13  modes  of  this  rectangle. 
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Tabic  3.  Results  for  the  Dominant  (Noncutoff)  Hj  Mode  of  the  Weakly  Guiding  2:1 
Rectangular  Guide.  Listed  values  of  p2  were  obtained  after  a  1  2-,  3-,  or  4- 

term  truncation  of  the  field  expansion.  The  field  Intensity  at  two  points  on  the 
perimeter  is  also  given.  The  field  intensity  at  p  0  Is  1 


d 

P 

2 

Intensity  at 
#  -  0  8  •  */2 

1  Term 

2  Terms 

3  Terms 

4  Terms 

3.0 

0. 9512 

0.9732 

0. 9726 

0.9711 

0.007 

0.023 

2.  5 

0.  9333 

0.  9628 

0.  9607 

0.  9597 

0.009 

0.031 

2.0 

0. 9035 

0.  9444 

0. 9406 

0.  9402 

0.  013 

0.047 

1.  5 

0.  8488 

0.  9074 

0. 9022 

0.  9021 

0.023 

0.075 

1.0 

0.  7354 

0.  8170 

0. 8116 

0. 8118 

0.047 

0.  142 

0.  5 

0.  4456 

0. 5149 

0. 5122 

0.  5124 

0.  155 

0.355 

Furthermore,  for  the  first  mode  of  each  family  (that  is.  R*),  the  single-term 

2  1 

truncation  provides  accuracies  in  I*  ranging  from  0.04  to  0.  14.  As  In  the  case 
of  the  square,  the  l -term  truncation  results  shown  in  Table  4  were  obtained  by 
uttlmng  only  the  dominant  term  in  the  appropriate  determinantal  Eq.  (.SO)  or  (51). 

41  Mbptiral  Crass  Swims 

The  same  calculations  which  were  performed  for  rectangular  guides  were 
repeated  using  circular  (ft  l)  and  elliptical  (ft  >  1)  cross  sections.  Although  the 

modes  for  the  generalised  ellipse  are  described  mathematically  by  the  same  four 

I  IV 

forms  listed  in  Eq.  (41),  they  are  labelled  Efl  through  E^  In  order  to  Identify  them 
as  pertaining  to  the  ellipse. 

The  propagation  constants  of  the  circular  guide,  which  arc  the  same  (that  is. 
degenerate)  for  the  sine  and  cosine  modes,  derive  from  single-term  truncation. 
They  are  shown  over  the  range  Os  d  <  2.  5  in  Figure  9.  As  ft  increases  from  1  to 
2,  the  sine  and  cosine  modes  split  and  shift.  The  results  for  ellipses  with  ft  •  1.2, 
1.  5  and  2  are  depleted  in  Figures  10,  It,  and  12,  respectively.  These  modes  are 
qualitatively  very  similar  to  those  of  their  rectangular  counterparts  shown  In 

V  V 

Figures  4.  5.  and  6.  In  fact.  If  one  makes  the  correspondence  E  —  R  .  then  the 

11  W 

order  in  which  the  first  doren  modes  of  the  elliptical  guides  reach  cutoff  la  nearly 
the  same  as  the  order  for  the  rectangular  guides.  Moreover,  the  actual  values 
of  d  at  cutoff  are  similar  for  the  two  shapes  for  all  ft.  Table  5  gives  the  actual 
values  of  d  at  cutoff  for  ellipses  with  l  9  ft  <  2  and  also,  for  comparison,  for  the 
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Figure  11.  Mode*  of  the  Weakly  Guiding  Elliptical  Guide  with  R  1.  S 


Figure  12.  Modes  of  the  Weakly  Guiding  Elliptical  Guide  with  ft  »  2 
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Table  5.  Values  of  ;l  at  Cutoff  for  Elliptical  Guides  with  Different  Aspect  Halloa 
and  for  a  Rectangular  Guul.'  with  ft  2.  Asterisks  indicate  where  the  designations 
KV  and  Ey  should  be  interchanged 


Values  of  (3  at  Cutof 

For  Ellipse 

For 

Rectangle 

ft  1 

ft  1.2 

ft  1.5 

ft  2.0 

ft  2.0 

0 

0 

0 

0 

0 

0.  78 

0.  68 

0.  58 

0.  48 

0.  42 

0.  78 

0.  74 

0.  72 

0.  67 

0.  60 

1. 23 

l.  12 

1.01 

0.  92 

0.  77 

1.23 

1.07 

0.  89 

0.  72 

0.  64 

1.  28 

1. 22 

1.  18 

1.  15 

l.  08 

1.  65 

1. 48 

1. 30 

1.  13 

0.  95 

l.  83 

l.  59» 

1.48* 

1.  37* 

1.  18 

1.  77 

1. 46* 

l.  22* 

0.  96* 

0.  87 

1. 77 

1.  72 

1. 68 

2.03 

1. 82 

1.  52 

1. 20 

1.  12 

2.  03 

1.  83 

1. 60 

1.35 

1.  14 

rectangular  guide  with  ft  2.  Furthermore,  the  correspondence  between  the 

^  ni 

and  R|'  modes  extends  to  the  field  configurations,  which  comprise  very  similar 
patterns. 

The  qualitative  similarities  between  modes  of  the  rectangular  and  elliptical 
guides  approach  numerical  agreement  as  the  aspect  ratio  becomes  large,  in  fact, 
for  H  10  the  first  few  modes  of  either  shape  guide  are  virtually  indistinguishable 
for  j?  >  1.  5.  In  general,  the  propagation  constants  for  either  differ  by  more  than 
0.  03  onlv  very  close  to  cutoff. 
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Figure  15.  Mode*  of  the  Weakly  Guiding  Cusp-Shaped  Guide  *tth  ft  1.  N  0.3 

There  are  many  atmtlaritiea  between  the  elliptical  and  rectangular  guldea,  so 
it  i*  no  aurpriae  that  all  the  meaaurable  chars*  -leristlcs  of  the  intermediate  super- 
elltpttcal  ahape  (the  promulgation  ronatanta,  field  conf igurat Iona,  cutoff  values, 
order  of  the  no-dos,  dependence  on  ft,  convergence  of  the  aeries)  fall  between 
those  of  the  two  former  shapes.  In  fact,  for  N  2  and  ft  l,  the  superelltpsc 
(Figure  13)  is  closer  to  the  square  than  the  circle  in  all  respects,  and  for  N  5 
(Figure  14)  it  is  practically  indistinguishable  from  the  square.  For  ft  1  the 
modes  all  shirt  in  a  sm*>oth  and  orderly  manner  as  N  changes  from  1  to  30. 

The  convergem  e  of  the  series  is  marginally  better  for  the  elliptical  shape 
than  for  the  superell ipses  and  rectangles. 

The  cusped  shape  (Figure  15)  is  quite  different  in  appearance,  and  the  charac¬ 
teristics  of  the  modes  reflect  this.  The  order  in  which  the  modes  cut  oft  is  quite 
different,  as  are  the  propagation  constants  themselves. 


One  can  demonstrate  that  the  individual  modes  "evolve"  as  the  aspect  ratio 
changes.  Figures  1«  to  20  show  the  propagation  constants  of  certain  selected 
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Kljfurc  18.  The  Hj'  Mode  of  the  Wraklv  Guiding  If  e<-t  angular  Guide  aa  R  C'hangea 
from  1  to  2 
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modes  versus  ,d  for  different  rectangular  croaa  sections,  ranging  front  a  square 

to  an  elongated  rectangle.  The  behavior  depicted  here  ia  typical  of  all  modea  In 

that  cutoff  occurs  at  smaller  $  as  ft  increases.  The  evolution  of  these  same  five 

nu>des  lor  a  continuous  change  in  ft  la  displayed  in  Figure  21  for  two  values  of  0. 

For  certain  modes  the  dependence  of  P*  on  ft  is  much  stronger  than  for  most 

others.  This  is  especially  evident  in  the  it!,  mode  (see  Figures  16  and  211  and  the 
tl  iv 

H.  (Figures  17  ami  21).  The  Rj  mode  (Figures  16  and  21),  which  Is  degenerate 

with  the  Hj'  for  ft  1,  shows  only  a  very  small  dependence  on  ft.  This  difference 

can  tie  understood  when  one  refers  to  the  field  diagrams  for  the  various  modes  as 

a  n 

shown  in  (Joell.  *  The  I).  field  has  two  maxima  along  the  x-axis  (the  long  direc- 

1  IV 

turn)  and  one  along  the  y-axla.  whereas  the  reverse  holds  true  for  R.  .  As  the 

II  1 

x -dimension  is  lengthened,  the  two  R  lobes  spread  a  considerable  distance  from 

1  IV 

the  center  along  the  x-axts,  whereas  the  Rj  maxima,  although  squeezed  Some¬ 
what,  do  not  shift  so  dramatically.  In  other  words,  the  greater  distortion  of  the 

It!1  field  results  in  a  greater  change  in  the  propagation  constant  as  ft  increases. 

*  I 

The  same  holds  true  for  R_  and  It.  (three  maxima  in  one  direction,  one  in  the 
II  iv  J 

other)  and  Rj  and  R-,  (four  and  one).  As  a  contrast.  Figures  1 F*.  20,  and  21 
display  examples  of  moderate  dependence  on  ft,  as  evidenced  by  the  l(|  and 
modes,  respectively.  In  F'lgure  21  one  can  see  much  more  vividly  such  effects 
as  the  splitting  of  the  degeneracy  in  the  rV*Rj'  mode  as  the  cross  section  departs 
from  a  square,  and.  the  crossing  of  the  It!,  propagation  constant  curve  with  those 
of  the  Hj"  and  Hj'  modes  as  ft  increases. 

One  can  also  see  from  figures  such  as  17,  19,  and  21  tliat  for  very  elongated 
rectangles  (ft  •  '>)  the  propagation  constant  is  indej«endent  of  ft  except  at  small  /J. 
This  suggests  that  the  greatest  field  strengths  occur  near  the  center  and  that  the 
intensity  at  the  ends  is  considerably  less.  This  is  indeed  true  except  near  cutoff 
(small  ,d)  where  the  field  spreads  throughout  the  guide  and  extends  into  the  sur¬ 
rounding  region. 

In  tracing  the  "evolution"  of  a  mode  as  ft  increases  one  can  avoid  confusion  by 
examining  the  field  configurations.  It  happens,  for  example,  that  the  R^'  mode 
for  ft  <  l.  2  corresponds  to  the  R1^  mode  for  ft  »  1.5,  and  vice  versa.  That  la, 
the  field  patterns  show  that  the  same  mode  which  derives  from  the  second  root  of 
the  equation  del  ( f J>  0  when  the  shape  is  close  to  square,  corresponds  to  the  third 
root  when  the  cross  section  becomes  elongated.  I*ut  another  way.  Individual  modes 
may  "cross  over  as  ft  changes  (just  as  they  may  aa  3  changes).  However,  the 
designation  R*  is  defined  for  a  particular  guide  (ft  fixed),  so  corresponding  modes 
can  be  labelled  differently  for  different  ft.  In  Table  5  a  similar  circumstance 
occurs  for  ellipses. 

For  the  first  mode  of  each  family  (namely  R^S  we  find  that  the  series  converges 

^ 

nearly  as  well  for  large  aspect  ratios  (R  10)  as  for  cases  with  R  close  to  unity. 
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and  four-term  truncation  provided  reasonable  accuracy  for  the  second  modes  us 

well.  However,  for  the  higher  modes  of  each  family,  5.  6,  and  even  7  terms 

2 

were  required  before  the  values  of  1’  would  he  obtained  to  within  0.03  for  such 
elongated  cross  sections. 

Figure  22  exhibits  the  same  information  for  the  ellipse  us  was  displayed  in 

Figure  21  for  the  rectangle.  As  before,  one  can  see  how  sensitive  each  mode  is 

11  IV 

to  change  in  aspect  ratio,  and.  how  a  mode  (for  example,  the  Fj  -Ej  )  which  is 
degenerate  for  a  circular  (A  1)  cross  section  splits  as  the  cross  section  becomes 
elliptical. 

1.0  r 


$  •  1.0 
£  •  to 


1.0  1.5  2j0 

ff 

2 

Figure  22.  Variation  of  P  with  Change  in  Aspect  Hatio  for  Some  of  the  First 
Few  Modes  of  the  Weakly  Guiding  Flltptic  al  Guide  at  Two  Values  of  (1 

Figures  Ifi  to  22  demonstrate  how  different  modes  are  affected  by  a  change 
In  the  aspect  ratio  (ftl  for  a  fixed  rectangular  or  elliptical  shape.  This  was 
accomplished  by  keeping  N  constant  in  Eq.  (fi0>  and  varying  ft.  Of  course,  sim¬ 
ilar  types  of  analyses  could  be  carried  out  for  changes  in  shape  (N>  aa  the  aspect 
ratio  remains  fixed. 
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It  is  * r tl  known  tliat  the  problem  of  finding  tlve  propagation  constants  of  a 

waveguide  having  a  particular  geometry  la  identical  to  the  problem  of  finding  the 

energy-  levela  of  a  quantum  mechanical  potential  well  whose  potential  haa  the  aame 

geometry  aa  the  waveguide.  The  Helmholtz  equation  tnalde  the  well,  given  by 

IT2  ♦  !  V I  -  Ik!),  O.  where  'v|  ka  the  well  depth,  'e|  la  the  energy  level  depth 

and  t  la  the  wave  function,  la  the  quantum  mechanical  counterpart  of  our  Eq.  (2). 

The  llalm holt z  equation  outaide  the  well,  given  by  *  |e|)»  •  0,  la  the  quantum 

counterpart  to  our  Kq.  (3).  Thua.  for  the  caaes  ronatdered  here,  there  ia  a  one* 

to -one  correapondence  between  the  appropriate  component  of  the  electromagnetic- 

field  <K,  B,  or  •)  and  the  quantum  wave  function  (*),  between  |v|  -  |e|  and 

yf  k*  -  k^.  between  !e!  and  y?  k^  -  ki!.  and  between  |v|  and  y^  ♦  y? 

llg  *  g  *  *  * 

k^  -  k2.  Moreover,  the  variation  in  the  propagation  conatanta  of  the  varloua  mode, 
with  changing  ft  but  fixed  ahape  (aee  Ftgurea  16  and  22)  la  auggeative  of  the  rhangca 
which  the  energy  levela  would  experience,  including  aplitting  of  degeneraclea,  aa 
a  perturbing  influence  la  applied  to  the  potential  well.  For  the  waveguide  caaea 
treated  in  Figuren  16  to  22,  the  role  of  the  perturbat ion  ia  played  by  the  changea 
in  the  gutde'a  aapect  ratio  aa  ita  ahape  remaina  fixed.  Alau,  aa  waa  mentioned 
above,  it  ia  poaaible  that  the  role  of  the  perturbation  could  be  played  by  the  changea 
in  the  gutde'a  ahape  for  a  conatant  aapect  ratio. 


(V  C — | - lilh  (Mm  *«i 

In  order  to  teat  the  efficacy  of  our  method  we  compare  our  reaulta  with  pre- 
vtoualy  publtahed  data. 

For  guidea  of  aquare  and  rectangular  croaa  aectiona,  there  are  mainly  the 
2 

reaulta  of  Goell.  His  method  conxiata  of  matching  field  componenta  across  the 

waveguide  boundaries  at  a  large  number  of  points  and  then,  similar  to  our  method, 

obtaining  a  determinants!  equation  whose  solutions  are  the  propagation  constants. 

From  these  aoluttona  he  then  determines  the  coefficients  in  hts  field  expansions. 

1 

Mar  cat  ill,  employing  a  similar  method,  analyzes  a  few  caaes  of  the  rectangular 

1 2 

guide  having  various  aapect  ratios.  F.yges,  on  the  other  hand,  utilizing  the  per¬ 
turbation  techniques  mentioned  in  Section  6,  Investigates  Some  of  the  lower  modes 
for  guides  having  small  values  of  ft.  Both  authors  report  data  which  are  in  general 
agreement  with  (Joel  I . 

Some  of  (Joell's  reaulta  from  the  first  six  modes  are  shown  in  Figure  3  for  the 

square  and  in  Figure  6  for  the  ft  •  2  rectangle.  They  were  taken  as  carefully  as 

possible  from  his  Figures  16  and  17  (square  and  ft  2  rectangle,  respectively). 

taking  note  of  the  fact  that  his  f)  differs  from  ours  by  a  factor  of  2.  Also  in 

Figure  «  we  show  the  results  from  the  perturbation  prescription  for  the  dominant 

<R,)  mode  of  the  ft  -  2  rectangle.  (It  turns  out  that  the  latter  method  predicts 
*  2 

values  of  P  for  the  dominant  mode  equal  to  those  of  a  circular  guide  of  the  same 


X  V 

amil.  For  the  dominant  mode  (Kj‘  j  in  Goell's  nutation)  there  la  no  discernible 
difference  between  hia  curve  and  our  a,  and  only  amall  differences  for  the  perturba¬ 
tion  analysis  results.  For  several  other  modes  there  are  nearly  uniform  differ- 

o 

ences  of  less  than  0.07  in  P  between  Goell's  results  and  ours.  That  is,  the 
curves  are  adjacent  and  parallel,  and  certainly  correspond  to  each  other.  For 
the  square  guide,  the  same  tyi>e  of  results  have  also  been  found. 

Many  qualitative  features  which  can  be  seen  in  our  Figures  3  and  <5  —  such  as 
the  splitting  of  It*,  and  lt*(  (square)  near  i  utoff  and  the  crossing  of  hJ'  and 
(rectangle)  —  are  also  present  in  Goell's  curves,  although  there  are  numerical 
discrepancies,  imly  with  the  It***  <E*'  ',)  mode  of  the  rectangle  are  there  dis¬ 
crepancies  greater  than  0.  1  in  P*. 

With  the  square  guide,  all  the  odd  modes  are  degenerate  In  pairs  (sine  and 
cosine)  because  of  symmetry.  (See  Figure  3  and  Table  2.)  Moreover,  except 
for  R*  and  H***,  those  even  modes  in  which  the  first  term  dominates  are  degener¬ 
ate  in  pairs  away  from  cutoff.  This  can  be  seen  clearly  in  Figure  3  for  the 

H_-R*  modes.  Furthermore,  if  the  .(-axis  were  extended  to  larger  values,  then 
*  J  in  hi 

this  degeneracy  would  be  more  ap|>arent  for  the  It.,  -lt.(  modes,  as  well  as  for 

higher  modes  which  are  not  shown. 

As  a  result  of  the  foregoing,  one  could  take  a  linear  combination  of  the  field 
expansions  for  two  of  these  modes  and  derive  an  intensity  pattern  cor  responding 
to  neither  of  the  individual  modes,  hut  rattier  to  the  combination.  In  fact,  with 
our  technique  this  would  automatically  occur  If  one  did  not  recognize  the  simplifi¬ 
cations  implied  by  the  symmetry  and  therefore  effect  the  separation  of  the  field 
expansions  at  the  beginning.  In  Goell's  Figure  5,  intensity  patterns  for  those 
modes  which  are  degenerate  reflect  this  combination  of  what  appear  —especially 
for  the  even  indices  —  to  be  distinct  modes.  For  rectangular  cross  sections 

(ft  »  l)  our  field  intensity  patterns  agree  with  Goell's  (see  his  Figure  6). 

1  17 

For  guides  of  rircular  cross  section,  there  are  the  works  of  (iloge  and  Yeh. 
All  of  our  results  for  this  cross  section  agree  very  well  with  theirs.  Yeh.  how¬ 
ever.  also  reports  data  for  guides  of  elliptical  cross  section.  A  compariscst  oi 
our  results  with  his  was  frustrated  by  what  may  be  either  notattonal  difficulties  or 
a  real  discrepancy.  The  problem  is  that  Yeh  lists  what  appear  to  be  two  different, 
but  degenerate,  modes  for  the  circular  guide;  ^HEjj  and  ^HE.j  In  his  notation. 
Then,  as  the  circle  deforms  into  an  ellipse,  he  shows  each  of  these  modes  becom¬ 
ing  a  nondegenerste  mode,  with  the  propagation  constants  of  one  mode  being 
greater  In  magnitude  than  those  of  the  circular  guide  and  those  of  the  other  mode 
being  leas.  However,  as  can  be  seen  In  Figure  22,  we  find  Just  one  dominant 
mode  (Ej)  for  the  circle  which  does  not  split  into  two  modes  as  ft  increases.  In 
addition,  our  results  for  this  single  noncutoff  mode  agree  very  we’l  with  a 

17.  Yeh.  C.  (197fi)()pt.  and  Quant.  Elect.  S:43. 
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12 

pei-turbation  calculation  based  on  Kyges.  Furthermore,  Yeh's  data  for  an 
elltpt leal  guide  having  an  aspect  ratio  of  2.  1«4  do  not  correspond  at  all  to  our 
results  for  either  an  ft  2  ellipse  (plotted  in  Figure  12)  or  an  ft  2.  104  ellipse. 
Hence,  are  are  unable  to  correlate  our  results  with  his. 

In  view  of  the  fact  that  our  results  for  the  rectangle  and  square  agree  quite 
closely  with  those  of  Cioell  and  the  fact  that  his  point  -matching  method  requires 
the  evaluation  of  a  determinant  of  order  4N.  with  3  <  N  <  9,  the  economy  of  our 
methtid  becomes  obvious.  Furthermore,  unlike  the  differential  equation  approach 
utilized  by  lloell,  our  integral  representation  technique  avoids  the  necessity  of 
having  to  use  a  coordinate  system  dictated  by  the  guide's  cross-sectional  shape 
as  well  as  the  s|>ecial  functions  characteristic  of  that  system.  This  means  that 
in  the  case  of  the  elliptic  guides,  for  example,  one  need  not  l>e  forced  to  use 
elliptic  cylinder  coordinates  nor  the  complicated  Mathleu  functions  associated 
with  such  coordinates. 

9  CONf.U  SION 


The  objectives  of  this  paper  have  hern  (II  to  present  a  general  method  by 
which  the  equations  for  the  propagation  constants  of  filler  waveguides  of  arbitrary 
cross  sectional  shapes  may  be  written  in  a  form  suitable  for  efficient  computation; 
(2)  to  specialize  the  method  for  weakly  guiding  fibers;  and  (3)  to  present  limited 
results  for  the  latter  case  to  demonstrate  the  correctness  and  efficiency  of  the 
approach.  The  extension  of  this  technique  to  coupled  waveguides  is  currently  in 
progress. 
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